Gaussian State Description of Squeezed Light 
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We present a Gaussian state description of squeezed light generated in an optical parametric 
oscillator. Using the Gaussian state description we describe the dynamics of the system conditioned 
on homodyne detection on the output field. Our theory shows that the output field is squeezed only 
if observed for long enough times or by a detector with finite bandwidth. As an application of the 
present approach we consider the use of finite bandwidth squeezed light together with a sample of 
spin-polarized atoms to estimate a magnetic field. 
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INTRODUCTION 



Quantum mechanical squeezing of optical fields repre- 
sents a means to improve precision measurements below 
the standard quantum noise limit for optical detection. 
These measurements exercise a significant back-action on 
the probed system, and to assess the achievement of a de- 
tection scheme we need a formalism that can deal with 
light-matter interaction and measurement induced state 
reduction continuously in time. Now, a continuous wave 
beam of light is described by infinitely many modes, for 
example in time or frequency domain, and the quantum 
state of the light field and of the system interacting with 
the beam is in general too complicated to be fully ac- 
counted for in a Schrodinger picture representation. In 
many quantum optical problems with constant or peri- 
odic driving Hamiltonians, it has been possible, however, 
to provide solutions in the Heisenberg picture for the rela- 
tionship between Fourier transformed frequency compo- 
nents of the field and system observables. Unfortunately 
this well-established technique does not apply in conjunc- 
tion with measurements on the joint system acting locally 
in time and hence affecting all frequency components of 
the observables at each measurement event. The present 
paper paper does not provide a solution to this general 
problem. Instead, we shall demonstrate that for a specific 
dynamics restricted to a specific class of states, the so- 
called Gaussian states, a significant reduction in the num- 
ber of parameters needed to fully characterize the system 
enables a complete description. The squeezed light pro- 
duced in an optical parametric oscillator (OPO) is in such 
a Gaussian state, implying that the field is fully char- 
acterized by the first-order and second-order correlation 
functions of the field variables. The interaction with an 
atomic system may destroy the Gaussian character, but 
we shall restrict our attention to optical interaction with 
a large collection of atoms through an effective collective 
atomic observable, which may in turn be well described 
by a Gaussian quantum state. The Gaussian state for- 
malism p. LJ y| was recently employed [J, la, 6] for the 
off-resonant Faraday rotation-like interaction between a 
continuous beam of light and an atomic ensemble. To 



describe the interaction with a continuous wave of light, 
we proposed 0. 0. Q to treat the beam as a sequence of 
short segments of light incident on the atoms. In the in- 
teraction, each light segment acquires some entanglement 
with the atomic sample and causes a modification of the 
atomic state when the light segment is probed after the 
interaction. The description of the incident optical beam 
is simple if the state of the field factorizes in components 
corresponding to each short segment of the beam. This is 
indeed the case for a coherent state of light, representing 
a normal laser beam. Realistic sources of squeezed light, 
on the other hand, have a finite bandwidth of squeezing 
which implies that correlations exist between the field ob- 
servables at different times. Here, we extend the Gaus- 
sian formalism to the case of the continuous output from 
an optical parametric oscillator (OPO). We reproduce 
the properties of the squeezed optical beam which are al- 
ready known from standard quantum optics treatments, 
and we apply our formalism to the example of atomic 
magnetometry. 

The paper is organized as follows. In Sec.[n] we recall 
some results of the standard treatment of squeezing in an 
OPO. In Sec. IIII1 we present the Gaussian state descrip- 
tion of this system. In Sec. lIVI we turn to magnetometry. 
The Larmor precession of an atomic sample caused by an 
unknown magnetic field is probed by optical Faraday ro- 
tation, and the value of the magnetic field is gradually 
determined as measurement data are accumulated. The 
use of a squeezed light source improves the magnetome- 
ter for probing times larger than the inverse bandwidth 
of squeezing. In Sec.El we conclude. 



II. GENERATION OF SQUEEZED LIGHT IN 
AN OPTICAL PARAMETRIC OSCILLATOR 



In this section, we present a simple model of squeezed 
light generation in a cavity with a non-linear medium 
which is pumped by a classical pump beam at frequency 
2w c , twice the cavity resonance frequency, and giving rise 
to creation and annihilation of pairs of photons by the 
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Hamiltonian 0] (we use h = 1 throughout) 

Hint = ig(a f2 - a 2 ) = g(x c Pc + PcX c ) (1) 

where a* and a are the creation and annihilation opera- 
tors for the light inside the cavity, and where the canon- 
ical conjugate variables are 



Xc = 7f( a + aT ) 



(2) 



We express the Hamiltonian in a frame rotating with the 
cavity resonance frequency ui c , and consider the dynamics 
in this rotating frame. 

In the absence of losses, the Heisenberg equations of 
motion 



x c (t) = 2gx c (t) 
p c (t) = -2gp c (t) 



(3) 



can be solved straightforwardly, leading to an exponential 
squeezing of the p c - variable and an accompanying anti- 
squeezing of the £ c -variable, which maintains a constant 
value of the uncertainty product. 

This model produces a squeezed state of a single light 
mode inside the cavity, and such states are subject to 
detailed analysis in most text books on quantum optics. 
Here, however, we aim at applications of squeezed light 
and hence we are interested in the squeezing properties 
of the light that leaks out of the cavity. This light prop- 
agates out of the cavity into a continuous beam, which 
corresponds to a continuum of modes in frequency space. 
We thus replace one of the perfectly reflecting cavity mir- 
rors with a mirror with a small transmittance, which will 
lead to a loss of the cavity field with rate T. The result- 
ing intra-cavity field state can be found in many different 
ways, but for our purpose it is sufficient to note that the 
cavity mirror acts as a beam splitter for the intra-cavity 
field and for the vacuum field (a; p h,in>Pph,in) incident on 
the cavity, see Fig. At the partly transmitting mirror 
the incident field is reflected into the output field. The 
output field is a linear combination of the reflected inci- 
dent field and the transmitted intra-cavity field. Imagine 
an incident beam segment of duration r, short enough 
that the intensity transmitted at the mirror and the field 
amplitude built up by the Hamiltonian can be treated 
to lowest order in r. We can then iterate the Heisenberg 
equations of motion for the intra-cavity field and the out- 
put field from the cavity and we obtain 



x c (t + t) = (£ + 2gr)x c (t) + VTTX phM (t) (4a) 

Pc (t + T) = (£ - 2 9 T ) Pc (t ) + VlVpph.in (t) (4b) 

£ph,out(^ + t) = -Vfrx c (t) + S,x ph:in (t) (4c) 
P P h,out(£ + r) = -Vfrp c (t) + £pph,m(i). (4d) 

where £ 2 = 1 — Tr denotes the probability for the seg- 
ment to be reflected by the mirror. This quantity is very 



close to unity, and consequently £ ~ 1 — Tt/2. The ex- 
pressions H4aH4d|l are of course equivalent to the ones ob- 
tained by the conventional input-output formalism 0,0], 
with the last terms in Il4al4l>ll having the characteristic 
properties of Wiener noise increments in the limit of small 
t. Since we assume that the input field is in the vacuum 
state, Eqs. H4al4hl) can be solved directly for the vari- 
ances of the intra-cavity field quadratures, starting from 
the vacuum at t = 0, and taking the r — > limit 



, , 1 T - Age^ r ' 4 9)t 

Var(;Cc)= 2 ~g 

ir + 4 ff e-( r + 4 ») t 
Var(Pc)= 2 TTTg ' 



(5) 
(6) 



If 4g < T, we see that these equations approach steady 
state for large times t. Since we are interested in oper- 
ating the OPO in a regime where steady state can be 
obtained, we assume from now on that 4g < T. The light 
inside the cavity is still squeezed as expected, but it is 
entangled with the emitted light, and hence it is not in a 
pure state and also not in a minimum uncertainty state. 

In the conventional input-output description, by a 
Fourier transformation to frequency space, the equa- 
tions Cfll become algebraic equations, and the output 
field operators in frequency space are expressed as linear 
combinations of the input operators at the same frequen- 
cies but with frequency dependent coefficients All 
moments of the field annihilation and creation operators 
have trivial expectation values in the vacuum state. If 
uj denotes the difference between the optical frequency 
and the cavity resonance frequency uj c , we have for ex- 
ample the following expression for the normal ordered 
expectation value of the output field when the system 
has reached steady state (remembering Ag < T) 



2r 5 



;6(W + U>'). (7) 



(i-2.g) 2 + ^ 



The Lorentzian frequency dependence implies a temporal 
correlation between the light emitted at different times, 
which is due to the common origin in the intra-cavity 
field. The field at a single instance of time is obtained by 
a Fourier transformation of the expressions in frequency 
space. This will involve all frequencies, also the ones far 
from the cavity resonance and hence outside the band- 
width of squeezing. Consequently, one will not observe 
squeezing properties if one observes a light field in a time 
interval shorter than ~ 1/T. Integrating the signal over 
a finite time interval T, corresponding to detection of the 
variable 
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t+T 



= —= / x(t')d£ 



(8) 
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yields a quantity with normal ordered expectation value 
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(9) 



If we use that (xj) — ('■ %t '■) + we see that for short 
times ({T — Ag)T <g; 1), the output field has the standard 
noise of vacuum, whereas integration over a longer time 
interval yields 



Var (:tt) 



l (r + 4g) 2 
2 (r - Agf 



(10) 



The corresponding variance for the pr component is ob- 
tained by replacing g by —g in the above expressions, i.e., 
in the long-time limit the emitted field is described by a 
minimum uncertainty state. 

The prediction of the noise properties of xt and pr 
should of course be in agreement with the ones observed 
if one carries out a homodyne measurement to detect 
these quantities, but it is important to remember that 
during such detection, the dynamics of the system will be 
different, and it is not clear how to modify the relations 
in frequency space between the intra-cavity and output 
fields as the detection takes place in real time. 

In the next section, we introduce the Gaussian state 
formalism which allows an effective real-time treatment 
of the production and probing of squeezed light. 



III. GAUSSIAN STATES 



A. General Formalism 




FIG. 1: Generation of squeezed light by an optical paramet- 
ric process pumped by a classical field at 2uj c . The figure 
shows three field segments in the vacuum state which enter 
the cavity, where a non-linear medium generates squeezing. 
Four segments of light are shown propagating away from the 
cavity. 

The linear transformation between the states of the 
cavity field and a segment of light initially incident on 
the cavity, and eventually propagating away from the 



cavity <@J is easy to deal with, because a state which is 
initially Gaussian in the field variables will remain Gaus- 
sian at later times. Taking an initially empty cavity and 
the incident vacuum field, will hence lead to Gaussian 
states at all later times. We want to consider the con- 
tinuous emission of light by the cavity, and we therefore 
imagine one segment of light after the other leaving the 
cavity, see Fig.nj and all field quadratures being given by 
a multi-mode Gaussian distribution. A Gaussian state 
is fully characterized by the mean value vector of all 
canonical variables, which we arrange in a column vec- 
tor y, with m = (y) and the covariance matrix 7 where 
jij = 2Re((y l - {yi)){y - (Vj))) • If the output field is dis- 
cretized in N segments m has dimension (27V + 2) with 
2N effective x and p variables for the output field and 
2 variables, x c , p c , for the cavity mode. The covariance 
matrix 7 has dimension (27V + 2) x (27V + 2). These fi- 
nite objects are of course far easier to deal with than the 
full TV + 1 tensor products of infinite dimensional Hilbert 
spaces. In practice the formalism can be made even sim- 
pler if we assume that the output beam is detected right 
after it is emitted from the cavity, and hence the quan- 
tum state of each light beam segment is destroyed and 
only the classical output value is retained, while the next 
segment emerges from the cavity. Let us consider the in- 
teraction between a single incident segment of light and 
the intra-cavity field, and let us write the linear transfor- 
mation of the four field variables y = (x c ,p c , a; p h,p p h) T 
as follows 



Sy 



(11) 



where the elements of the 4x4 matrix S follow directly 
from the transformation Under this transformation, 
the mean value vector m and the covariance matrix 7 
transform as 



m(t + r) = Sm(t) 
7(i + r) = S 7 (t)S T . 

We write the 4x4 covariance matrix as 



7 



C 7 B 7 , 



(12) 
(13) 



(14) 



where A 7 is the covariance matrix for the intra-cavity 
field variables, B 7 is the covariance matrix for the beam 
segment in the continuous beam, and C 7 represents their 
mutual correlations. An advantage of the Gaussian de- 
scription is that the back-action on the residual system 
due to measurement may be accounted for explicitly. If 
we measure the variable x p h, due to their mutual cor- 
relation, we learn something about the intra-cavity x c 
variable, i.e., its variance decreases, and simultaneously, 
to fulfill Heisenberg's uncertainty relation, Var(p c ) in- 
creases. Following Refs. Q, 0, HJ, we have the explicit 
update formula for the intra-cavity field covariance ma- 
trix after homodyne detection on the beam segment 



(7rB 7 7r) C 



7' 



(15) 
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where ir = diag(l, 0), and ()~ denotes the Moore-Penrose 
pseudoinverse. This result does not depend on the ac- 
tual outcome of the measurement. The latter however, 
affects the mean values of the intra-cavity field variables. 
The beam segment has disappeared from the treatment, 
but to treat the interaction with the next segment we 
build the covariance matrix (fill) describing the intra- 
cavity field and this new segment with 



B, 



1 
0. 



(16) 
(17) 



corresponding to an incident vacuum state with no cor- 
relation with the cavity field yet. We propagate the sys- 
tem according to Eq. ll.H . and we implement the effect 
of the subsequent measurement by Eq. J15II . The contin- 
uous production and probing of the beam is obtained by 
repetition of the above steps, and we may, in the limit 
of small time increments, derive a differential equation 
for the intra-cavity field A 7 . This differential equation is 
of the general non-linear matrix Riccati form (see, e.g., 
Ref. [ll| and references therein) 

A T (i) = G — DA 7 (<) - A 7 (i)E - A 7 (t)FA 7 (t), (18) 

where the 2x2 matrices G, D, E, F are all derived from 
the expressions l(T3j) and ltl5|l . As shown in Ref. 0] the 
solution for A 7 can be expressed in terms of the solutions 



of two coupled linear matrix equations: 
where W = -DW + GU and U = FW 



f EU. 



wu 



B. Squeezing properties of the intra-cavity field 

Applying the above general formalism to the squeezed 
light problem, we derive the Riccati equation itTsl) for 
the intra-cavity field covariance matrix, conditioned on 
the homodyne detection of the output field, and find the 
following matrices 



G = 


(0 0) 

\ o r ) 


F = 


(TO) 
\ o ) 


D = 


( -2g-V/2 
V 


E = 


D. 





2g+r/2 



(19) 



Without detection, the beam and the intra-cavity field 
are entangled, and, as we noted above, the intra-cavity 
field state, regarded as a trace over the unobserved emit- 
ted field degrees of freedom is a mixed state and not a 
minimum uncertainty state. If we perform homodyne de- 
tection on the emitted field we find the same variance of 
above, but the variance of p c changes to the value 



1 



Var(p c ) = - 



r-4 5 



2r-4 5 e-( r - 4 s)*' 



(20) 



In this case Var(x c ) • Var(p c ) = 1/4 and we have a min- 
imum uncertainty state of the intra-cavity field at all 



times. If 4<7 < T, we reach steady state for large times, 
and the variances then read 



Var(a; c ) 
Var(p c ) 



1 r 

1 r - 4. 9 

2 r ' 



(21) 
(22) 



In Fig. we show how the variance of the Gaussian 
variables inside the cavity depends on time both with 
and without measurements on the output beam. 




T 

2 4 6 8 

Time t in units of F _1 

FIG. 2: Variances of the cavity variables x c and p c as a func- 
tion of time. We use T = 2tt x 6 x 10 6 s' 1 and g = 0.2T 
which are realistic experimental parameters for OPO's [T^ |. 
The variances of x c with and without homodyne detection of 
the x p h variable of the output field are identical and shown by 
the upper dashed curve. The full and the dotted curves show 
the variances of p c without and with homodyne detection of 
the output field, respectively. 



C. Squeezing properties of the emitted beam 

1. Collective observable for many light segments 

We now turn to the squeezing of the output beam. 
As discussed in Sec. |nj there is no squeezing if we only 
consider small time intervals. To study the correlations 
between different individual segments we define the fol- 
lowing collective operators 

JY 



lT " VlvE x P h i ' 

i=l 
N 



(23) 



i=l 



where T = tN is the accumulated time in N segments 
each of duration r and where the field variables of the i th 
segment are retained in the formalism. In appendix EI 
we calculate the variances of these quantities. The result 



5 



is given in Eq. IIA13JI and reads 
1 



(V - Ag)(T + AgYT 



2T(r - 4 9 ) 3 
-32Tg + 32Tge { - r/2+29)T 



(24) 



The result for Var(pT) is obtained by replacing g with 
—g. If we let T — > 0, we obtain Var(p T ) = | showing 
that there is no squeezing if we only consider short times. 
If, on the other hand, we let T — > oo, we obtain 



Var(xr) 



l (r + 4g) 2 
2 (r - Agf 



(25) 



These results are in full agreement with the ones ob- 
tained by the usual quantum optics treatment discussed 
in Sec. El (see Eq. P)) l. 



2. Finite Bandwidth Detection 

An alternative way to extract the squeezed component 
of the emitted beam, is to use a frequency filter, that 
selects the frequency range of interest. The modelling of 
such a detector involves a second cavity, in which the light 
segments enter and the intra-cavity field in the second 
cavity builds up. The squeezed beam contains photons in 
the relevant frequency band, but not only intensity builds 
up in the detecting cavity, we also expect the intra-cavity 
field to show squeezing properties. The variables used in 



2u) c , 





FIG. 3: Proposed setup for the characterization of the spec- 
trum of squeezed light from an OPO (to the left). The squeez- 
ing properties of the single- mode field accumulated in the fre- 
quency tunable cavity to the right are determined (see text). 

a Gaussian treatment of this problem, corresponding to 
the two cavity fields and the propagating beam segment, 
are y = (x Cl ,p Cl , x C2 ,p C2 , a; ph ,p ph ) and the Heisenberg 
equations of motion are obtained by a simple extension of 
the expressions used already in the case of a single cavity, 
where we replace T with T±. The second cavity is used 
to model the finite bandwidth detection. It has a decay 
constant r 2 and a tunable cavity resonance frequency 
lo c + 6. In our frame rotating at to c the field variables in 
the second cavity obey the equations 



Eqs. H2ril27JI the field incident on the second cavity is the 
output field from the first cavity, cf. Fig. Due to the 
physical separation L of the two cavities and the finite 
speed of light the field variables in Eqs. H4cl4d|l should 
in fact have been delayed by L/c, but since we are ad- 
dressing the steady state properties of the system we can 
solve Eqs. Il4cl4dl) with the same time arguments. The 
output field from the second cavity is described by equa- 
tions similar to Eqs. I|4cl4dll . but they will not be needed 
in the following. The detuning 5 of the second cavity can 
be scanned, and the squeezing parameter of the intra- 
cavity variables x C2 , p C2 reflect the spectral properties of 
the output beam from the first cavity. 

Fig. shows the eigenvalues V m j n and V max of the 2x2 
covariance matrix for the probing cavity as fu nction of 
the detuning with respect to oj c . In panel 4(a) the prob- 
ing cavity has a damping rate r 2 comparable with the 
one of the OPO cavity, i.e., the intra-cavity field builds 
up with a memory time shorter than t he tim e needed to 
see the full effect of squeezing. In panel |4(b)| we use a de- 
tector system with narrow bandwidth, the cavity builds 
up light over a longer time inte rval, and the degree of 
squeezing is clearly larger than in 4(a) The Riccati equa- 



tion can be solved analytically, and for 5 = we obtain 



1 (IV 

2 (ri 



4g) 2 + (Ti - 4 5 )r 2 



4 5 ) (Ti 



4 5 ) 



(28) 



Here V m i„ is obtained from V maJC by replacing g with —g. 
The insert shows V m \ n and V max as function of IV For 
large r 2 , the second cavity is equally fed by a wide range 
of frequency components, and the variance is dominated 
by the vacuum uncertainty: Vmin = K max = 1/2. If T2 = 
we obtain 



i(r 1 + 4< 7 ) 5 



2 (Ti - 4. 9 )2 

1 (£1 - 4g) 2 
2(r 1+ 4 5 )2 



(29) 
(30) 



x C2 (f + r) = (1 - T 2 t/2)x C2 (t) + iSrp c 



p C2 (t + t) = (1 - r 2 r/ 2)p C2 (t) - iSrx c 



r2Tir phi0ut (t) 

_ (26) 

^2TPph,out(t) 

(27) 



which equal the long-time integrated amplitudes QSJ. 

We note that the calculations here were significantly 
easier than in the case where we treated a large number of 
light segments simultaneously. This is because the mode 
of the second cavity in practice integrates the incident 
field over time and stores the contribution of many short 
beam segments in a single set of variables. We believe 
that this is a useful model of realistic finite bandwidth 
detectors, and that the approach can be used quite gener- 
ally to investigate how finite optical bandwidth detection 
affects the sensitivity of metrology and the entanglement 
and spin squeezing of atomic samples. 



IV. MAGNETOMETRY WITH SQUEEZED 
LIGHT 



where a; p h,out, Pph.out are the quadrature variables for The purpose of introducing the Gaussian state formal- 
the field leaving the first cavity, cf., Eqs. Ij4cl4d|l . In ism is to provide a theoretical approach, that allows a 
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Detuning 5 in units of Ti 
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-1 -0.5 0.5 1 



Detuning S in units of I?i 

(b) 

FIG. 4: The variances V m i n and V mlK of the field inside the 
probing cavity as function of the detuning 5 of this cavity 
with respect to lo c in units of the decay width Fi of the OPO 
cavity. We have used r x = 2tt x 6 x 10 6 s _1 and g = 0.2Fi 
as in Fig. 2. In (a) T 2 = Ti and in (b) T 2 = Ti/25. The 
lower and upper parts of the insert in (b) show Vmin and 
Knax respectively as functions of the bandwidth of detection 
r 2 for 5 = 0. 



treatment of the interaction between light and an atomic 
sample in the regime where the quantum state of the 
atoms changes both because of the interaction itself and 
because the continuous measurements of the light field 
after the interaction teaches the observer about the state 
of the atoms. This measurement induced back-action on 
the quantum state of atoms plays a role in atomic mag- 
netometry > it was used to spin squeeze atomic gasses 
and to entangle pairs of gasses |"15l |. and it recently 
played an important role in the realization of an atomic 
memory for light [ifj- Since probing with squeezed light 
potentially is more precise, it was proposed in Ref. jj| 
that magnetometry would also benefit from the use of 
squeezed light, and a simple model with ultra-broad band 
squeezing indeed suggests improvement by precisely the 
squeezing factor on the B-field uncertainty. We will now 



use magnetometry as an example to show how we can ef- 
fectively treat the probing of atomic systems with a real 
squeezed optical field with finite bandwidth. 

It is possible to estimate a magnetic field by a polariza- 
tion rotation measurement of an off-resonant light beam 
passing through a trapped cloud of spin- 1/2 atoms, see 
Fig. El All the atoms are assumed to be polarized with 
their spin along the x direction. We assume that the 
B-field component of interest is directed along the y di- 
rection, and hence it causes a Larmor rotation of the 
atomic spin toward the z axis. This in turn leads to 
a mean magnetization of the sample along the z direc- 
tion, which will cause a Faraday rotation of the linear 
polarization of an optical field propagating through the 
sample. As the polarization rotation is proportional to 
the atomic spin component, and this is proportional to 
the B-field, the B-field estimated by the measurement is 
trivially obtained. We wish to address the error bar, i.e., 
the standard deviation on our estimate of the field as a 
function of the measurement record. The Gaussian state 
description which operates explicitly with the variances 
and covariance elements of the physical quantities is ideal 
for this analysis. 

The gas of trapped spin- 1/2 atoms is described by a 
collective spin operator J = (1/2) ^""i where <Xj are 
the Pauli spin matrices. The atoms are initially pumped 
such that they are polarized along the x axis and J x can 
be treated as a classical variable (J x ) = N at /2 where N at 
is the large number of atoms. The two other projections 
of the spin, J y and J z obey the commutation relation 
[J y , J z ] = iJ x which may be rewritten as [af a t)Pat] = i 
for th e effe ctive position and momentum variables x at = 
Jy/V( J x)' -Pat = Jz/y/(Jx)- The uncertainty is easily 
shown to be minimal in the initial state and, hence, the 
state pertaining to x at and p at is Gaussian. 

Note that the OPO cavity produces a squeezed vacuum 
state; if this field is linearly polarized along the y-axis, 
it may be mixed on a polarizing beam splitter with a 
classical x polarized field to yield the field appropriate 
for polarization rotation measurements. 
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FIG. 5: Setup for estimating a B-field. In the cavity, we 
generate squeezed light which is linearly polarized along the 
z axis. We mix this field at an asymmetric beam splitter with 
a strong x polarized beam. The light then passes through a 
gas of polarized atoms, causing a rotation of the field polar- 
ization towards the z axis. 

The light beam propagates along the y axis and is 
linearly polarized along x such that its Stokes operator 
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(S x ) = N p h/2 = 3>t/2 is classical with N p ^ the number 
of photons in a given segment of the beam and $ is the 
photon flux. The two remaining Stokes vector compo- 
nents corresponding to the difference in photon numbers 
with linear polarization along directions at 45 and 135 de- 
grees with respect to the x axis, and with left and right 
circular polarizations, respectively, have vanishing mean 
values, and they satisfy a commutator relation similar to 
the collective atomic s pin. Accordingly, f or the effective 
variables x ph = S y /y/(S x ), p ph = S z /y/{S x ), we have 
[£ph,P P h] = i for iVph S> 1 and the initial coherent state 
of the field is a minimum uncertainty Gaussian state in 
these variables. 

The effective Hamiltonian for this part of the system 

is 



Ht = Ky/rp at p ph + [irByXaf 



(31) 



The characteristic atom-light coupling is k = 

AAct "v^i^x)^ where d is the atomic dipole mo- 
ment, lj is the photon energy, A is the detuning of the 
light from atomic resonance, and A is the area of the 
the light field. T he co upling between the B-field and the 
atoms is /x = /3y/ (J x ) where (3 is the magnetic moment. 

The classical B-field that we wish to estimate is treated 
as a random variable with a broad Gaussian probabil- 
ity distribution. Hence both the B-field, the atomic 
cloud, and the incident light pulse of duration r are 
Gaussian variables. We arrange these in the vector 
y = {By, x aU Pat, x P h, P P h) T ■ The Larmor precession in 
time r and the interaction between the atomic sample 
and the beam segment then leads to the linear transfor- 
mation ifTTJl of the variables 0, 0] given by 



\ 





1 Ky/T 

-,ut 10 

Ky/? 1 

V o o o o i / 



(32) 



The initial covariance matrix is 70 = 
diag[2 Var(i? ), 1, 1, 1, 1]- After application of the 
matrix S, and the operations l(T3jl and l(T5)) representing 
the polarization detection of the optical field, corre- 
sponding to a homodyne detection of the variable x p h, 
the atomic variables and the B-field become correlated, 
a new beam segment enters in the standard coherent 
state as described by ltl6|l and ifTTJt . and the evolution 
proceeds. In the limit of short beam segments, the evo- 
lution can be replaced by a Riccati differential equation 
for the 3x3 covariance matrix for the atoms and the 
B-field, and this equation can be solved analytically Q. 
The variance of the B-field is 



Var(B(*)) 



Var(5o)(K 2 f + 1) 



which yields precisely the error on the estimate of the B- 
field. We note that the uncertainty of the field strength 
decreases as 1/ (N at t 3 ^ 2 ) and not as l/y/N at t as one might 
expect from standard counting statistics arguments. This 
improved precision is due to the squeezing of the atomic 
spin during the probing process. 



A. Squeezed light 

In Refs. Q we modelled the use of squeezed light by 
introducing the squeezing parameter r such that Eq. 11 fill 
is replaced by 



i^V Var(B )i 4 
6 



f«V Var(B )i 3 + n 2 t 



t^oo K 2 /I 2 t 3 



1/r 
r 



(34) 



i.e., every beam segment enters the interaction in a 
squeezed state. Going though the calculations we find 
that k 2 should be replaced with K 2 r in Eq. l(33jl : the 
B-field estimate is improved. 

As noted in Ref. this treatment of a squeezed beam, 
in the limit of small r, is only valid if the squeezing band- 
width is infinite. The squeezing properties of the beam 
from an OPO, however, only reveal themselves if a nar- 
row frequency component is selected, or if the field is in- 
tegrated over times longer than the inverse bandwidth of 
squeezing, which are certainly longer than the infinitesi- 
mal t employed in the continuous limit, where the Riccati 
equation is solved. 

The full probing may well take longer than the inverse 
bandwidth, and one would hence expect that one still 
benefits from the squeezing in this longer time limit. We 
shall verify this assumption by a calculation in which we 
treat the probing with the field coming out of our OPO 
cavity in the full Gaussian formalism. 

The example serves as a model for how to consider 
other atomic probing schemes with realistic squeezed 
light sources. We treat as Gaussian variables the B-field, 
the atomic variables, the intra-cavity field, and a single 
segment of light y = [B z , x at ,p at , x c ,p c , x ph ,p ph ) T . The 
beam segment enters on the cavity mirror in the vac- 
uum state, it is reflected off the mirror with some squeez- 
ing and some entanglement with the partly transmitted 
intra-cavity field, it interacts with the atoms, and finally 
it is detected by homodyne detection, causing a moderate 
change of the joint covariance matrix for the B-field, the 
atoms, and the intra-cavity field. The transformation to 
lowest order in r of the variables is now given by 



S = 



/I 














\ 




1 




















1 kVTt 





-k^Jt 











£ + 2 3 t 





Vtr 














£ - 2.9T 













-VTv 










\o 


— Ki/t 





-Vrv 








(33) 



(35) 
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with £ = 1 — IY/2 as introduced in Eq. J3J. Note that 
this matrix combines the elements present in the trans- 
formation of the field components alone 10} and the B- 
field-atom and light-atom interaction H32I1 . Again the 
beam segment is inserted in its vacuum state 111 (ill , and it 
is probed by homodyne detection leading to the update 
formula l(T5|l . The bandwidth is taken care of by the 
intra-cavity field which establishes the necessary correla- 
tion between beam segments detected at different times. 
In the continuous limit we find the corresponding Riccati 
equation, and its solution provides the variance of the 
B-field as a function of time as shown in Fig. 03 The fig- 
ure shows both the results without squeezing, with finite 
bandwidth squeezing, and the simple infinite bandwidth 
result (|33l) with a simple squeezing parameter r applied 

to each segment. We take the value r = [r^if^ 3 > corre- 
sponding to the long-time limit of Eq. il25l) , and we see 
a good agreement for long times between the two curves 
for squeezed states. We also see, that the finite band- 
width curve is an improvement with respect to the case 
of non-squeezed light, but that we have to probe for a cer- 
tain time on the order of the squeezing bandwidth before 
we see the effect of squeezing. Indeed, the finite band- 
width curve is to a good approximation simply delayed by 
-4g){r+4g) 2 compared with the infinite broad-band 
squeezed light curve. 

The analytical result for Var(B) is very lengthy. For 
small times t we get the result without squeezing as can 
be seen in Fig. 03 and for large t the result is exactly 
the same as in the infinite bandwidth case Vax(B(t)) = 
^■i^ K -i t i if we identify the squeezing parameter by r = 

(r+4 3 ) 2 
(r-4 ff ) 2 • 




20 40 60 

Time t in units of 



FIG. 6: Variance of the B field as a function of time. We use 
the same value of g and T as in Fig-EJ and k 2 = 1.83 x 10 6 s _1 
and fx = 8.79 x 10 4 (spT) _1 . The dashed line is without 
squeezing, the full line is with squeezed light generated in 
a cavity, and the dotted line is with the squeezing parameter 
r. 



V. CONCLUSION AND OUTLOOK 

In summary, we have presented a Gaussian state de- 
scription of the light from an optical parametric oscillator 
and its interaction with large atomic samples. The treat- 
ment is very effective, because the state of the parts of the 
beam that have just left the OPO cavity can be treated 
as a single mode, corresponding to a short beam segment, 
and after the interaction, the segment can be eliminated 
from the formalism. Here, we presented the dynamics 
when the field is probed by homodyne detection, and it 
is turned into classical information; if the beam propa- 
gates away without detection, it may be traced out of 
the formalism, which is an even simpler operation in the 
Gaussian formalism, since the corresponding rows and 
columns in the covariance matrix should just be removed. 
Finite bandwidth effects are included in the treatment 
by retaining the quantum state of the intra-cavity field, 
which is also a single field mode, i.e., at the price of 
adding a single pair of canonically conjugate variables 
(x c , p c ), which in the Gaussian formalism is done by 
adding two extra rows and columns to the covariance 
matrix. 

The use of squeezed light holds the potential to im- 
prove spin squeezing, entanglement, and precision prob- 
ing, and we demonstrated such an improvement in the 
case of magnetometry compared with the infinite band- 
width case, we also showed how the finite bandwidth of 
squeezing manifests itself as a time lag before the im- 
provement is obtained in agreement with the observation 
that squeezing is only present in a light beam, if one in- 
tegrates a sufficiently long part of the beam. 

The method described is fully general, and further 
studies can be carried out along the same lines on other 
proposals involving squeezed light. It is readily gener- 
alized to incorporate more atomic systems, more field 
modes, non-degenerate OPOs, and as we showed also fi- 
nite detection bandwidth can be modelled by the addi- 
tion of auxiliary modes. Finite bandwidth of the light 
sources and of the detection system may also play non- 
trivial roles in conjunction with decay and decoherence 
which set an upper limit to the degree of entanglement 
obtained in gasses [l7j . 

Finally we note that squeezed light has been proposed 
as an ingredient in various quantum information proto- 
cols, such as teleportation pj|, as a source of heralded 
single photons [lil l2fT | , as a resource in continuous vari- 
able quantum computing and error correction [2J. In 
many of these protocols, an elementary analysis is given 
in terms of single mode fields, where indeed, a full time 
and frequency dependent analysis would be more appro- 
priate. 
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APPENDIX A: SQUEEZING OF OUTPUT FIELD 

To calculate the variance of xt and px defined in 
Eq. l|23jl we use the Gaussian description. 




PIG. 7: The figure shows how we label the beam segments 
which are accumulated in xt- We only obtain squeezing if 
we observe many beam segments, not if we only observe one 
segment. 

The initial variables are yi = (x 

c i Pa x ph 1 j Pph x ) and 
the initial intra-cavity field and incident vacuum segment 
covariance matrix is 



7i = 



/On 0,12 ON 

ai2 a 2 2 

10 

\ 1/ 



(Al) 



The associated transformation matrix is given by Eq. <@J 



2gr 
£ - 2gr 

-VTv 
o -vTv 



r- 






Tt 


C ) 



(A2) 



After the 



where £ = 1 — IY/2 as introduced in Eq. 
interaction ■y 1 = S^jSf . We now build the dynamics 
recursively by inserting two rows and columns between 
the second and third row and column in 7. In this way 
we represent the subsequent incident vacuum segments 

by 



7jh 



{7fc}(l:2,l:2) 


v {7fc}(3:2/c+2,l:2) 





1-2x2 





{7fe}( 

{7fc}(3: 2/c+2,3:2fc+2) . 



r(l:2,3:2fc+2) 





The transformation matrix is 
S fc 



Si 

1 







(2fc-2)x(2fe-2) 



(A3) 



(A4) 



and 7 fc = S k j k S k . Eqs. IJA4IA3I) are now in- 
serted, and the number of variables grows with time 
as we get more and more light segments, yjy = 

{%c,Pc, x ph N j Pph N j ■ ■ ■ > x ph 1 i Pph ± ) 

If Q42 = then every second element in 7 is zero and 
7 can be rewritten on block diagonal form with similar 
x p h and p p h blocks. The system of equations for the a; p h, 
Yn = {x Cl Xph 1 , ■ • ■ , ^phjy ), variables can be written as 

'l - IY/2 + 2gr vTV 
-a/IV l-IY/2 
1 




'A k CV 
Ik = I 1 
>c£ B fc/ 



(A6) 



(A7) 



where A k is a real number, C k is a 1 x (fc— 1) row vector, 
and Bfe is a (fc — 1) x (fc — 1) matrix. 

From this we find the recurrence equations 



A k+1 = (1 - rY/2 + 2gr) 2 A k + Tt 



C T - 
^fe+i — 



(A8) 



rr(l - IY/2 + 2gr)A k + \/Tt(1 - IY/2) 



(1 - Tr/2 + 2g T )C{ 



B 



fe+i — 



fTrA k + (1 - Tr/2) 



-v/ivcr 



IYC fe 
B fc 



(A9) 
(A10) 



which can be solved, and the variance of xt is found to 
be 



^ N N 

Var(x T ) = — Cav(xi,Xj) 
i=i j=i 

Trl- a 2JV 



= an 



2N 1 



+ 



rva 4 



1 - a 



JV 



1-a 



2^ 



2V(1 - tt)\l-a 1-a 2 
rr 2 Tr 2 l-a 2W 



2(1-0*) 2A^(l-a 2 ) 2 
+ - Tr/2) 2 -1Y(1- IY/2) 



1 



1 - a 



IY(1 - IY/2) 1 - a 



N 



+ 



N (1-a) 2 

Tr 2 a / l-a N 

N 



N(l-a 2 )(l-a) 



1 — a 



1 — a 



2N 



N 



I — a 2 



1-a 
1 — a 



(All) 

where a = 1 - Tr/2 + 2gr. If we let T = Nt and then 
let t — > then 



Vax^ x _ {(r-4 ff)( r + 4^ 



r 



2T(r-4 5 )3 
-4T(T + 8.g) + 4o u r(T - ig) 

-sr[ ail (r - 4 5 ) - (r + 4. 9 )] e (- r / 2 + 2 ^ 
-4T[r-on(r-4 5 )] 



,(-r+4 S ) 



1- 



(A12) 



In the T — > 00 limit, the first term dominates, and the 
expression for Var^r) does not depend upon an. 
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In steady state we may insert an = r _ 4g 
and we obtain 



m Eq. for the variances of the quadrature components. 



Var(xT) 



and 



Var(p T ) 



1 



(r-4g)(r + 4 5 ) 2 T 



2T(r-4g)3 
-32Tg + 32Tge^ r/2+2g)T 



(A13) 



(r + 4.g)(r-4 5 ) 2 T 



2T(r + 4 5 ) 3 
-32Tg - 32.g(r - Ag)e^ r l 2+2 3) T (A14) 

-64 5 2 e-( r+4 ^ T l 
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